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We investigated macroscopic quantum tunneling (MQT) of Bi2Sr2CaCu2 0j, intrinsic Josephson 
junctions (IJJs) for two device structures. One is a small mesa, which is a few nanometers thick with 
only two or three IJJs, and the other is a stack of a few hundred IJJs in a narrow bridge structure. 
The experimental results regarding the switching current distribution for the first switch from the 
zero-voltage state were in good agreement with the conventional theory for a single Josephson 
junction, indicating that the crossover temperature from thermal activation to the MQT regime 
for the former device structure was similar to that for the latter device structure. Together with 
the observation of multiphoton transitions between quantized energy levels in the MQT regime, 
these results strongly suggest that the observed MQT behavior is intrinsic to a single IJJ in high-Tc 
cuprates and is independent of the device structure. The switching current distribution for the 
second switch from the first resistive state, which was carefully distinguished from the first switch, 
was also compared with respect to the two device structures. In spite of the differences between the 
heat transfer environments, the second switch exhibited a similar temperature-independent behavior 
for both devices up to a much higher temperature than the crossover temperature for the first switch. 
We argue that this cannot be explained in terms of self-heating caused by dissipative currents after 
the first switch. As possible candidates for this phenomenon, the MQT process for the second 
switch and the effective increase of the electronic temperature due to the quasiparticle injection are 
discussed. 



I. INTRODUCTION 

Macroscopic quantum tunneling (MQT) in a current- 
biased Josephson junction is a key phenomenon exhibit- 
ing quantum mechanical behavior in systems of macro- 
scopic degrees of freedomjii^ii^ MQT is observed only at 
sufficiently low temperatures where thermal activation 
(TA) over the potential barrier is suppressed, while TA 
behavior is exhibited at temperatures higher than the 
relevant crossover temperature Tcr. Recently, the obser- 
vation of MQT has attracted great interest as a fun- 
damental measurement for manipulating the so-called 
"phase qubit" , which is regarded as one of the promising 
candidates for future quantum information processing.— 
In such applications, some conventional superconductors 
(typically, Al and Nb) have been used for the fabrication 
of Josephson junctions. Thus, most experiments regard- 
ing the observation of MQT and the operation of phase 
qubits were performed only at sufficiently low tempera- 
tures (- 0.1 K). 

The discovery of high temperature cuprate supercon- 
ductors has drastically changed the situation due to the 
high critical temperature Tc and the high critical cur- 
rent Ic of these superconductors, which are one or two 



orders of magnitude higher than those of conventional 
superconductors. Particularly, intrinsic Josephson junc- 
tions (IJJs) in Bi2Sr2CaCu20y (BSCCO) single crystals 
have attracted great interest since a crystalline stack of 
the superconducting Cu02 double layer and the insulat- 
ing Bi202 (or SrO) layer forms a perfectly smooth tunnel 
junction on an atomic scalei^ Indeed, recent experiments 
using a small stack of IJJs in a BSCCO single crystal have 
clearly indicated that MQT behavior can be observed up 
to nearly 1 Ki^i 

However, the use of stacked IJJs with almost identical 
Ic gives rise to a serious problem where it is impossi- 
ble to ensure that a switch to the first resistive state al- 
ways occurs at the same junction when the bias current 
is repeatedly increased to obtain the probability distri- 
bution of the switching current. In addition, the anoma- 
lous enhancement of the MQT rate for "uniform-stack 
switching" , where several tens of junctions in IJJs switch 
to the resistive states at almost the same time, remains 
unresolved)^ Very recently, to avoid this complication of 
using stacked IJJs, a surface IJJ formed on the top of a 
small mesa with five IJJs was proposed by Li et al,^ since 
Ic of the surface was significantly smaller than that of in- 
ner IJJs in the stack. The results from their experiments 
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regarding the surface IJJ, which indicate that MQT be- 
havior is observable only below ~ 0.1 K, seem to point to 
another problem: the properties of the surface IJJ might 
be strongly influenced by a thick normal-metal electrode 
which is in direct contact with the surface IJJ. Thus, pre- 
vious studies^'^'^ of MQT in stacked IJ Js and surface IJ Js 
have suggested that experiments involving the MQT phe- 
nomenon in IJJs should be performed more carefully by 
considering various influences originating from the device 
structure and the fabrication method of IJJs. 

In this paper, we present a comparative study of MQT 
in the following two types of IJJs. One is a small mesa, 
a few nanometers thick with only two or three IJJs, and 
the other one is a stack of a few hundred IJJs in a nar- 
row bridge structure. The probability distribution of the 
switch from the zero-voltage state to the first resistive 
state in the multiple-branched current-voltage character- 
istics was carefully measured for both types of IJJs, con- 
firming that the same junction switched to the resistive 
state in each event. We demonstrate that the MQT be- 
havior for the surface IJJ on the former mesa structure 
survives up to around 0.4 K, which is similar to Td- for 
the stacked IJJs in the latter bridge structure, in con- 
trast to the results obtained in the study performed by 
Li et al.^ We also observed microwave-induced multipho- 
ton transitions between quantized energy levels in the 
first switch for both types of IJJs, which were expected 
to appear in the MQT regime. The observation of both 
MQT and energy-level quantization (ELQ) in the first 
switch strongly suggests that these phenomena are in- 
trinsic to a single IJJ of BSCCO and are independent of 
the device structure. 

Moreover, we performed similar switching-current dis- 
tribution measurements for the second switch from the 
first resistive state to the second resistive state in the 
current- voltage characteristics. In spite of the differences 
in the heat transfer environments, the second switch for 
both devices showed a temperature-independent behav- 
ior, which was similar to MQT in the first switch, up 
to much higher temperatures (approximately an order of 
magnitude higher) than Tcr for the first switch. We dis- 
cuss this unusual behavior from various points of view. 
In contrast to the recent paper by Kashiwaya et al.^ the 
experimental results suggest that this behavior cannot be 
explained in terms of self-heating caused by dissipative 
currents after the first switch, as will be discussed below. 

This paper is organized as follows. In Sec.|TTl we briefly 
describe the conventional model of the escape from the 
zero-voltage state, which is based on the resistively and 
capacitively shunted junction (RCSJ) model. In Sec. Illli 
the preparation of both types of IJJs and the method 
employed in the switching-current distribution measure- 
ments are described. Section lTVl contains all experimental 
results and a discussion. Finally, we conclude this work 
in Sec. El 



II. MODEL 

The current- voltage characteristics of a current-biased 
single Josephson junction are successfully explained 
on the basis of the so-called RCSJ model, which 
was independently studied by Stewart and McCumbcr, 
respectively!^ In this model, the phase difference, 6, 
across the Josephson junction is described by the follow- 
ing equation of motion fii 

(fe I dO ^ I 

-dP^Q-,iit^''^'-Tr'^ 

where Ic is a critical current and Q^^ is a damping pa- 
rameter given by Qo = t^poRC, respectively, R and C 
are the shunt resistance and the shunt capacitance, re- 
spectively, and LOpo is the so-called Josephson plasma fre- 
quency given by 

-PO = (2) 

As is well known, the phase dynamics described by Eq.([T]) 
is equivalent to a damped motion of a particle with a 
mass of (?i/2e)^C in a tilted washboard potential, U{9), 
which is given by 

U{9) = -Ej{cose + je), (3) 

where E,j{— hlc/2e) is the Josephson coupling energy, 
and 7 = I / Ic- 

In the underdamped region (Qo > l)i where the 
current-voltage characteristics become hysteretic, the 
zero- voltage state corresponds to the trapped particle at 
the bottom of a potential well (that is, 9 is fixed at a 
certain value). When the washboard potential is tilted 
with increasing the bias current, the particle escapes from 
the well and starts rolling downhill (that is, 9 changes at 
a rate d9/dt — 2eV/h). The bias current at which the 
junction switches from the zero- voltage state to the finite 
voltage state is referred to as switching current, Isw 

In the absence of thermal and quantum fluctuations, 
such a switch always occurs at 7 = 1. However, when 
thermal (or quantum) fluctuations cannot be neglected, 
the switch is induced by a TA process over the poten- 
tial barrier (or an MQT process through the potential 
barrier) for a smaller bias current 7 < 1. The potential 
barrier is given hy^ 

AU{j) =i;,/[7(2sin"^7-7r) + 2cos(sin~^7)]. (4) 

The escape rates induced by TA and MQT processes are 
given by 
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respectively"^, where ujp is the current-dependent Joseph- 
son plasma frequency, 



LOp = Wpo(l - 7 ) 



2U/4 



(7) 



at(= 4/[^l + {QkBT/1.8AU) + 1]^) is a thermal prefac- 
tor, where Q {=Qo{l — j'^Y^^) is the current-dependent 
damping factor. Eq.® describes the escape rate induced 
by MQT without dissipation. At high temperatures, the 
TA process is dominant, while a crossover to the MQT 
process occurs at the crossover temperature, which is ap- 
proximately given by 
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Experimentally, as a function of the bias current, the es- 
cape rate, t(/)~^, is related to the probability distribu- 
tion, P{I), of the switching current through the following 
equatiop^^^ 



t{I)-' ^ In 
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P{i) 
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where A/ is a small current interval, the probability for 
the switch between / and / -f A/ is given by P{I)AI, 
and dl /dt is the ramp rate of the bias current. Thus, it 
is expected that the standard deviation of which 
is given by a{= ((/gw - (^w))^)^^^), is proportional to 
r2/3 in the TA regime since AC/(7) « £;j4v/2/3(l - 
7)'^/^ for 7 close to 1, while a is independent of T in 
the MQT regimci^ The fitting results of P{I) using Ic 
and an effective temperature at the escape event, Tosc, 
also suggest that Tosc agrees with the bath temperature, 
Tbath, hi the TA regime, while Tcsc becomes independent 
of Tbath in the MQT regime.?^ 

When microwaves are irradiated onto the junction in 
the MQT regime, the distribution of P{I) is divided into 
multiple peaks, which is attributed to the quantized en- 
ergy levels in the potential well. The spacing between the 
n-th and the n + 1-st energy levels, En,n+i, is obtained 
by solving the Schrodinger equation under the boundary 
condition that the eigenfunction becomes zero when its 
energy is equal to the potential energy for a free-running 
side of the barrier. Numerical calculations suggest that 
Eqi is roughly given by fiojp as a function of At/(7)i^ 



III. EXPERIMENTAL METHOD 
A. Preparation of IJJ devices 

In this work, we studied two types of IJJs in BSCCO 
single crystals, which were grown by the floating zone 
(FZ) or the traveling-solvent floating zone (TSFZ) 
methodi^ii^ and annealed to be optimal doped. 

First, a small mesa, a few nanometers thick, with a 
few IJJs referred to as "mesa-type IJJs" below, was fab- 
ricated as follows. A piece of the single crystal with a 



typical size of 0.5x0.5x0.01 mm"^ was glued to a sapphire 
substrate and cleaved in a vacuum chamber before evap- 
orating a silver film to reduce the contact resistanceii^ A 
small mesa structure with a height of a few nanometers 
was fabricated by using electron-beam (EB) lithography 
and Ar-ion milling techniques. The number of included 
junctions was confirmed by the observation of the I — V 
characteristics, as shown below. Note that IJJs of this 
type are usually measured with a three-terminal config- 
uration, as shown in Fig. [T{a). 

On the other hand, a few hundreds of stacked IJJs 
in a narrow bridge structure, which are referred to 
as "bridge- type IJJs" below, was fabricated as follows. 
First, four gold electrodes were deposited on another 
piece of the BSCCO single crystal with dimensions of 
1x0.1x0.01 mm"^. Next, a narrow bridge ~ 5 fim long 
and ^ 1 /im wide was fabricated between two voltage ter- 
minals by using the focused ion beam (FIB) technique. 
Then, the sample was rotated by 90 degrees in order to 
etch the side wall of the bridge. Finally, two slits were 
fabricated at a distance of ~ 1 fxm by using FIB etch- 
ing in order to separate a small stack of IJJs from the 
bridge structure. A typical scanning ion microscope im- 
age is shown in Fig. 1(b) The number of included junc- 
tions was estimated to be about 200 from the height of 
the stacks. Note that the bridge-type IJJs are measured 
with the usual four-terminal configuration, in contrast to 
the mesa-type IJJs. The details regarding the sample 
parameters are summarized in Table ID 



B. Measurement technique and apparatus 

The probability distribution of the switching current 
from the zero-voltage state (or the first resistive state) 
was measured down to 0.4 K by using a ^He cryostat, 
as followsi^ First, a bias current is ramped linearly with 
time by using an arbitrary-waveform generator (Agilent 
33220A) until the junction switches from the zero- voltage 
state to the first (or second) resistive state, and is reset to 
zero in order to wait for the next sweep. Next, as shown 
in Fig. 1(c), the ramped current and the voltage across 
the junction, which are amplified with a battery-powered 
preamplifier (Stanford Research Systems, SR560) or an 
operation amplifier (Analog Devices, AD524), are mea- 
sured as a function of time by using a digital oscilloscope 
(Tektronics, TDS5054) or a data acquisition board (Na- 
tional Instruments, NI6251). The switching current, Igw, 
is recorded when the voltage is beyond a threshold value, 
which is typically set to 5 % of the switching voltage in 
the first (or the second) switch above the zero-voltage 
state (or the first resistive state). Finally, the probabil- 
ity distribution of Igw is obtained by repeating a large 
number of measurements of Igw (typically 5000~ 10000 
events) . 

High-frequency noise through the four feed-lines con- 
nected to the IJJ sample is eliminated by using the fol- 
lowing three stage low-pass filters. The first stage is a 
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TABLE I: Various parameters of the fabricated IJJs. The lateral size and the number of IJJs for the mesa-type samples were 
determined by the mask size for EB lithography and the I-V measurements, respectively, while those for the bridge-type samples 
were determined by the image of the scanning ion microscope. Tc for each sample was determined by the zero-resistance, which 
was measured under a bias current far below the switching currents. Ic (or 1^"'^) was determined by the least square fitting of 
log(r~'^(7)) for the first (or the second) switch (See also Ref. [18]). Ta was estimated from Eq. (|8}, and the uncertainty in Tcr 
was about 30% due to the range variation of Up, er, and 7c. 
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Type 


size (/im^) 


N 


Tc (K) 


Ic (M) 




i!"" (M) 


M 


mesa 


2x2 


2 


74 


10 


0.4 


41 


N 


mesa 


2x2 


3 


43 


2.7 


0.2 




R 


mesa 


2x2 


3 


50 


6.3 


0.3 




S23 


bridge 


1.1x0.7 


~200 


70 


4 


0.6 




S301 


bridge 


1x1 


~200 


84 


21.3 


0.9 


23.8 


S302 


bridge 


1x1 


~200 


84 


12.1 


0.8 


22.3 




100/im 1//m 



(c) DitirAequisition Board 




FIG. 1: (Color online), (a) (top) Schematic representation of 
the bias lead configuration and (bottom) a micrograph image 
of the mesa-type IJJ device, (b) (top) Schematic configura- 
tion and (bottom) a scanning ion micrograph image of the 
bridge-type IJJ device, (c) Schematic representation of the 
measurement setup. 



temperature of decompressed '^He. The second stage fil- 
ter is a combination of a 1 Ml chip resistor inserted at 
1.6 K and the lossy coaxial cable, which was found to play 
a role similar to a RC low-pass filter. The final stage filter 
is a home-made powder filter consisting of a 1.8 m long 
twisted pair of Manganin wires, which is coiled up in a 
50 mm long copper tube filled with 100 mesh size SUS304 
powder. The room-temperature measurement performed 
by using a network analyzer (HP8753D) confirmed that 
the attenuation of this powder filter was more than 50 dB 
at frequencies above 1 GHz. 

The electromagnetic radiation from the environment is 
also carefully shielded as follows. First, the IJJ sample is 
mounted in a copper cell which is thermally anchored on 
the cold finger of a ^He pot. Next, both a /i- metal and 
a permalloy shield are prepared in order to reduce stray 
magnetic fields. Moreover, all experimental components 
which do not need an AC power supply, such as the load 
resistor for monitoring the current, the battery-powered 
preamplifiers, the three stage filters, the sample cell, and 
the liquid-helium cryostat, are placed inside a shielded 
room (Japan Shield Co. Ltd., ARll), while the other 
instruments, which need an AC power supply, are placed 
outside the shielded room. 

Another coaxial transmission line consisting of several 
0.085" coaxial cables and V-connectors was fed into the 
sample cell in order to observe the microwave irradiation 
effect up to 67 GHz. We also used a short superconduct- 
ing NbTi coaxial cable for the thermal isolation of the 
sample cell and a coaxial attenuator (10 dB or 20 dB) 
for the reduction of thermal noise in the transmission 
line. 



IV. RESULTS AND DISCUSSION 



low-pass TT-type feedthrough filter with a cutoff frequency 
of 1.9 MHz, which is operated at room temperature. We 
use a lossy stainless-steel coaxial cable (0.3 mm in diam- 
eter and 1 m in length) with the attenuation of 24 dB / m 
at frequencies above 1 GHz as the feed-lines to intervene 
between room temperature and 1.6 K, which is the base 



A. Properties of the first switch 

Figures[2fa) and[Hb) show the typical I—V character- 
istics of the mesa-type IJJs (sample M) and the bridge- 
type IJJs (sample S302), respectively. They were mea- 
sured at a temperature of approximately 5 K. The nega- 
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tive voltage offset associated with the measurement elec- 
tronics was subtracted in Figs. [IJa) andjU^b). As shown 
in Fig. [2lja), only two resistive branches were observed in 
the I — V characteristics for the mesa-type IJJs (sample 
M). No more branches were observed up to higher bias 
(which is not shown here), which clearly indicates that 
the total number of junctions included in this sample was 
only two. We also note that the voltage jump (~10 mV) 
for the first switch in this sample is still larger than the 
previously reported value for similar mesa- type IJJs.^ 

On the other hand, we found that the voltage jumps for 
the bridge-type IJJs were not always uniform, as shown 
in Fig. [ll^b). In the case of sample S302, the voltage 
spacing between the zero- voltage branch and the first re- 
sistive branch was slightly larger than the spacing be- 
tween the second and third (or the third and forth) resis- 
tive branches. Since the latter voltage spacing (roughly 
25 mV) is in good agreement with the previously reported 
values, ^'-S we consider that the former larger voltage spac- 
ing is attributed to the underdoping effect due to the 
residual Ga ions in the FIB etching process. A similar 
effect at the surface IJJs in the the mesa structure has 
been discussed by Zhao et ali^ In contrast to the first 
voltage jump, the next voltage spacing between the first 
and second resistive branches was too narrow to be ex- 
plained by the local overdoping. Rather, it is very similar 
to the behavior in a resistively shunted junction, implying 
that the junction which becomes resistive at the second 
switchings is locally shunted by a small resistance. 

We also found that the retrapping current (roughly, 
3 fiA) in the second branch was slightly larger than that 
in other branches. This feature is apparently similar to 
the behavior which was recently observed for the mesa- 
type IJJs with a small number of junctions (typically, less 
than 10)3ii2, Although You et al}^ suggested that this 
feature was attributed to the sheet critical current in a 
topmost layer of the mesa structure, we found that our 
result could not be explained by this scenario, since the 
observed difference of the retrapping currents between 
the first and the second branches was at least one or- 
der of magnitude smaller than the estimate of the sheet 
critical current (We roughly estimated it as 70 ^A using 
the value given by You et al.). Rather, the larger retrap- 
ping current in the second resistive branch is qualitatively 
explained by the small shunted resistance, since the re- 
trapping current is inversely proportional to the damping 
parameter, Qo{= UpoRC). Thus, we consider that the 
unusual features observed for the second resistive branch 
are attributed to the small shunted resistance. It is also 
important to note that these features are not common 
to the other bridge-type samples. Although the origin 
of such a small shunted resistance is still unknown, it is 
clear that the small shunted resistance does not affect the 
switching current measurements. 

As shown in Figs. [UJa) and[2l^b), it was observed for 
both samples that the switching current for the first (or 
the second) switch was definitely smaller than that for 
the second (or the third) switch. The large difference 
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FIG. 2: I ~V characteristics of (a) sample M at 5 K and (b) 
sample S302 at 4.5 K. 



between the first and the second (or the second and 
the third) switching currents allowed us to distinguish 
the corresponding switched junction from the rest of the 
junctions. We consider that the origin of the reduced 
switching current for the first switch is different in the 
mesa-type and the bridge-type samples. Regarding the 
mesa-type samples, the first switch is considered to occur 
at a surface junction, as reported by Li et a/4 However, 
their interpretation that the surface layer was intention- 
ally underdoped seems to be difficult to explain the ob- 
served small voltage jump. Rather, we consider that a 
kind of proximity effect plays an important role in reduc- 
ing the switching current at the surface junction.^" On 
the other hand, in the bridge-type sample, the residual 
Ga ions seem to reduce the switching current in the first 
and the second switch. Although the microscopic details 
of the influence exerted by such residual Ga ions are un- 
clear, a plausible explanation is that the residual Ga ions 
make nearby junctions slightly underdoped, leading to 
the reduction of the critical currents. 

These results clearly suggest that the quality of the 
fabricated IJJ devices is sufficiently high for investigat- 
ing the switching properties in the TA and MQT regimes. 
Based on these observations, the respective switching 
current distributions, P{I), for the first switch in both 
samples (M and S302) were measured for temperatures 
down to 0.4 K, as shown in Figs. [31(a) and[3l^b). We found 
that the experimental data (solid symbols) for P{I) be- 
low 4 K fitted the theoretical calculations (solid lines), 
where we used Ic and Tosc as fitting parameters for the 
single junction model in the TA regime and assumed sim- 
ply that at was 1.— Figures [31(c) and[Hd) show the tem- 
perature dependence of the standard deviation of P{I), 
cr(T), for several samples. The experimental results (solid 
symbols) of cr(T) were also compared with the calcula- 
tions (solid lines) in the TA regime. We observed that 
the behaviors of cr(T) for both mesa- and bridge-type 
IJJs were divided into the following three temperature 
regions; (i) the higher temperature region (above approx- 
imately 4 K), where (t{T) increased as the temperature 
decreased, (ii) the intermediate temperature region (be- 
tween approximately 4 K and 1 K), where cr{T) decreased 
as the temperature decreased, (iii) the lower temperature 
region (below approximately 1 K), where (t{T) began to 
saturate as the temperature decreased further. 
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Temperature ( K ) Temperature ( K ) 

FIG. 3: (Color online), (a) Switching current distribution, 
P(/), for the first switch of sample M. (b) Similar plots of 
P{I) for the first switch of sample S302. (c) Solid symbols 
are the temperature dependence of the standard deviations, 
(t(T), of P{I) for the mesa-type IJJs (sample M and N), while 
solid lines represent the calculation in the TA regime, (d) 
Similar plots of cr(r) for the bridge- type IJJs (sample S302 
and S23). 



As shown in Fig.[3ljc), we found not only that (t{T) for 
the mesa-type IJJs (sample M and N) was proportional 
to T^/-^ in the region (ii), but also that <j{T) was in good 
agreement with the calculation in the TA regime, where 
it is assumed that T^sc = T as well as that at = 1. This 
strongly suggests that the first switch for the mesa-type 
IJJs are successfully described by the conventional sin- 
gle junction model. The crossover temperature, T^^, was 
estimated from Eq.dS]) as 0.4 K and 0.2 K for sample M 
and N, respectively, as listed in Table [D Wc found that 
the estimated T^.^ was slightly lower than the tempera- 
ture below which criT) showed a saturated behavior in 
the crossover region from the TA process to the MQT 
process. Unfortunately, the reason for this small discrep- 
ancy between the two temperatures remains unknown at 
present. However, our results clearly suggest that the 
crossover from the TA process to the MQT process oc- 
curs in region (iii). Furthermore, the observed behav- 
ior in region (i) is qualitatively explained by the phase 
retrapping-diffusion processi^ 

The results for a{T) for the bridge- type IJJs (sam- 
ple S302 and S23) were found to be almost the same 
as those for the mesa-type IJJs, as shown in Fig. ^d). 
However, a small discrepancy in cr{T) between the exper- 
imental data and the calculations in the TA regime was 
observed in region (ii). At present, the following three ex- 
planations can be considered for this discrepancy. One is 
the contribution of the thermal prefactor, at, which was 
neglected in the above analyses. Since at is dependent 
on Q^, new fitting analyses including the contribution 
of at suggest that Q of the bridge-type IJJs was more 
than 250, while that of the mesa-type IJJs was 30 ~ 70. 
As discussed below, Q was estimated to be 100 for the 



bridge- type IJJs from the resonant peak under microwave 
irradiation. This suggests that the above fitting analyses 
tend to overestimate the Q value. The second one is the 
effect of large Josephson junctions, which becomes im- 
portant when the junction size is larger than the Joseph- 
son penetration depth, Aj. A previous study, in which 
a large mesa of IJJs (from 10x10 to 40x40 /im^) was 
utilized, suggested that the switching distribution was 
further broadened by the large junction effects.^"* By us- 
ing the fitting results in the TA regime and the assumed 
value of Aab(~ 250 nm), Aj was estimated to be about 
0.4 fim for sample S302, which is indeed smaller than the 
junction size (1 /im). However, we found that the esti- 
mated Aj(~ 0.8 /im) for sample M was also smaller than 
the junction size (2 /im), suggesting that this effect is un- 
able to explain the small discrepancy in cr{T) observed 
only in the bridge-type IJJs. As the third candidate, we 
found that the discrepancy in cr(T) became more promi- 
nent when the switching current for the first switch was 
closer to that for the second switch. This implies that 
the coupling effect between adjacent layers, which is ne- 
glected in the conventional single junction model, might 
be closely related to the observed discrepancy in cr(T). 

In order to confirm the realization of the MQT process 
more clearly, we measured P{I) under microwave irra- 
diation at the lowest temperature. In the MQT regime, 
where the energy levels in the potential well are quan- 
tized, the transition from the ground state to the first ex- 
cited state can be induced by microwave radiation whose 
frequency is equal to the energy level spacing in the po- 
tential well. The numerical calculation of the Schrodinger 
equation suggests that the energy level spacing, £'oi, 
between the ground state and the first excited state is 
roughly given by hujp, as described in the previous sec- 
tion. Since hujp/2TT easily exceeds 100 GHz in the case 
of IJJs, it is difficult to observe the transition across Eqi 
due to the single photon absorption process. However, 
the multiphoton process based on using a microwave fre- 
quency equal to one for integer Eqi is expected to work 
well even in the present case.— 

Figure Ufa) shows the measurements of P{I) for sam- 
ple S301 under microwave radiation with a frequency 
of 45.7 GHz. We observed that a resonant peak de- 
veloped in the lower current region of P{I) as the mi- 
crowave power increased, which corresponded to the 
MQT process from the excited state. We estimated the 
current-dependent plasma frequency '^p{= 2itx 136 GHz) 
at 7=19.0 /iA by using the values of the critical current 
/c=21.3 /iA and the junction capacitance C=39.9 fF, 
which is estimated from the junction size and the rel- 
ative permittivity er(~5)-^- Note that the estimated Wp 
is three times as large as the applied microwave frequency 
at 45.7 GHz, suggesting that a three-photon process was 
detected. 

Figure Sfb) shows an equi-intensity contour in the 
switching current versus the microwave power plane 
for three-photon absorption, which was measured at 
45.7 GHz for sample S301. Here, the microwave power 
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FIG. 4: (Color online), (a) P{I) for the first switch of sample 
S301 under several different levels of microwave irradiation 
(b) Equi-intensity contour in the switching current ver- 
sus microwave power plane for three-photon absorption, mea- 
sured at 45.7 GHz for sample S301. (cl)(c2) l/r(/) under 
microwave irradiation for samples S301 and R, respectively, 
(d) Plots of the applied microwave radiation versus the nor- 
malized resonance current for several samples (S301, S302, 
and R). Solid, dashed, dotted and dash-dotted lines represent 
the calculations for single-photon, two-photon, three-photon 
and four-photon absorptions, respectively, which are based on 
Eq. 0. 



is measured at the output port of the microwave synthe- 
sized sweeper. Since the attenuation of the microwave 
power in the transmission line, including some connectors 
and an attenuator, depends strongly on the microwave 
frequency, it is difficult to compare the power at the 
sample position between the different frequencies. We 
also plot the escape rate, r(/)~^, which was obtained 
from P{I) for sample S301 and sample R, as shown in 
Figs. HJcl) and Ill;^c2), respectively. Note that t(/)~^ 
for sample S301 corresponds to a three-photon process 
(a;p/27r=45.7x3 GHz), while that for sample R corre- 
sponds to a single-photon process (wp/27r=55 GHz). We 
found that the range of t{I)~^ for sample S301 slightly 
shifted toward the lower current region as the microwave 
power increased, in addition to the growth of the res- 
onant peak due to the MQT process from the excited 
state. This is in sharp contrast to the results for t(/)~^ 
for sample R, where no shift in the range of t(/)^^ was 
observed as the microwave power increased. We found 
that these behaviors were related to the fact that the mi- 
crowave AC current affects the switching characteristics 
of the junctions. In this case, the switching current 
distribution is strongly affected by the microwave power. 

In order to take into account the influence of the ac 
current on the switching current, we corrected the ap- 
parent difference between the switching current under 
microwave irradiation and that without microwave ir- 
radiation. As shown in Fig. Hid), we plotted the res- 
onant frequency, fp{= ujp/2Tr), which was observed for 



several multiphoton processes as a function of the cor- 
rected bias current at which At(/)^^/t(/, 0)~^[= 
-fshiftjPmw)^^ - t{I,0)^^}/t{I,0)^''-] showed a resonant 
peak. Here, /shift and Pmw indicate the apparent shift 
in the switching current which compensates the contri- 
bution of the ac current and the irradiated microwave 
power, respectively. The solid lines in Fig.HJd) represent 
the calculation based on Eq. ([7]), which is in good agree- 
ment with our experimental results. We also estimated 
the Q value for such resonance peaks by fitting them 
with a Lorentzian, which suggested that Q was 80^^110 
for sample S301 and 50^90 for sample S302, respectively. 
Regarding sample R, for lower rci(~0.3 K), we observed 
only a broad step-like feature in At{I)^^ /t{I,0)^^ for 
sample R even at the lowest temperature, in contrast 
to the results for samples S301 and S302 with higher Ta- 
Similar results have already been reported by Martinis et 
al^ suggesting that the measured temperature was not 
far below Tcr. 

It has been highlighted recently that the microwave- 
induced multi-peaked switching distribution and the 
Rabi-type oscillation can be explained by using the clas- 
sical model of Josephson junctionsi^ This scenario is 
likely to become more important for multilevel systems 
with larger dissipation. In order to distinguish between 
a two-level system and a multilevel system, it is useful 
to evaluate AU/hujp.^'^ In our samples, AU/hujp was es- 
timated to be about 5 for samples S301 and S302 and 
about 10 for sample R at the lowest temperature. The 
values for the former samples, S301 and S302, are typi- 
cal when observing Rabi oscillation in a two-level system, 
suggesting that our results shown in Figs.HJa) and|Hb) 
can be regarded as indicating quantized energy levels. 
On the other hand, the value for sample R is closer to 
that for a multilevel system. This behavior can be ex- 
plained by considering the fact that AU/htUp decreased 
as the bias current increased. Since the switching current 
decreased as the temperature increased near and above 
Tcr, the larger value of AU/Uuip for sample R suggests 
that the measured temperature is not far below Tcr, as 
discussed above. Thus, a clearer indication of ELQ for 
sample R is expected to be observed if the temperature is 
decreased further. Based on these results, we concluded 
that the indication of MQT and ELQ in the first switch 
was essentially common to both device types, which sug- 
gests that it was intrinsic to a single IJJ of BSCCO and 
independent of the device structure. 



B. Properties of the second switch 

We also measured the temperature dependence of P{I) 
for the second switch from the first resistive state in 
the same way. Figure [S] shows the plots of Tesc versus 
the bath temperature, Tbath, where Tesc was estimated 
from Eq. (O . We found that the decrease of Tesc for the 
first switch saturated below a certain temperature, T^^^ 
(roughly '--^l K), while that for the second switch satu- 
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FIG. 5: (Color online). Plots of the effective escape tem- 
perature Tesc versus the bath temperature Tbath for both the 
first and the second switch of (a) the mesa-type IJJs and (b) 
the bridge-type IJJs, respectively. The solid lines represent 
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FIG. 6: (Color online), (a) l/r(/) for the second switch of 
sample S302 at 0.42 K for different values of dl/dt. Corre- 
sponding pulse widths are 33.3 msec {dl /dt = 1 mA/sec) ~ 
0.333 msec {dl/dt = 100 mA/sec). (b) Similar plots of l/r(7) 
for different waiting times between current pulses. Solid, dot- 
ted and dashed lines represent the fitting to Eq.([5} with the 
current-independent rosc(=10.2 K), with Tesc oc / and with 
Tosc cx: respectively. 



rated below a much higher temperature, T^'^'^ (roughly 
~8 K). Surprisingly, these behaviors were found to be 
common to both types of IJJs. It should be noted that 
T^^"'^ was also much higher than the estimated crossover 
temperature, T^,^""^ (~0.8 K for the mesa-type samples 
and ^1.3 K for the bridge- type samples), for the second 
switch. Note that the saturation at higher temperature 
is not attributed to the heating of the contact resistance, 
which was less than 100 17, since its influence was not 
observed in the first switch. 

Similar behaviors have recently been investigated by 
Kashiwaya et al^ They argued that this unusual behavior 
was attributable to the self-heating effect in the resistive 
state after the first switch. However, we found that this 
interpretation could not explain our results successfully 
since the self-heating effect was strongly dependent on 
the heat transfer environment in the device structure, 
which was significantly different in the mesa-type IJJs 
and the bridge-type IJJs. The problem of self-heating 
for the mesa structure has been numerically investigated 
by Krasnov et al.^^ They estimated the magnitude of the 
self-heating by using the two-dimensional solution of a 
heat diffusion equation for a small circular mesa on top 
of a semi-infinite single crystal. They also suggested that 
the contribution from the metallic electrode on top of the 
mesa and that from an insulating layer separating the 
electrode from the bulk crystal played an important role 
as parallel heat channels. On the other hand, it seems 
that the heat diffusion process in the bridge-type IJJs 
with a few hundred junctions is quite different from that 
in the mesa-type IJJsi^ For instance, the heat diffusion 
process along the stacked direction (that is, along the c- 
axis of the crystal), becomes more important since the 
stack of IJJs is far away from the remaining part of the 
bulk crystal. In addition, there is no effective parallel 
heat channel since the stack of IJJs is surrounded by 
vacuum and is located at a sufficient distance from the 
electrode. This strongly suggests that the magnitude of 



the self-heating in the bridge-type IJJs should be much 
greater than that in the mesa-type IJJs. 

If we assumed that the thermal resistance, Rth, was in- 
dependent of the dissipated power, Pdis, we can estimate 
the temperature rise, Arh(= ^th x ^tiis), at a current 
of the second switch using the recent experiments of the 
self-heating which were measured at a larger dissipated 
powerr-^^^i^i For instance, Verreet et alS^ estimated that 
i?th ~160 K/mW at Pdis ^150 /iW for the mesa-type 
IJJs. This suggests that ATh for sample M is ^--^0.14 K for 
Pdis =0.9 /iW in our measurements at the lowest temper- 
ature. On the other hand, Wang et alM- estimated that 
Pth ~450 K/mW at Pdis ~70 for the bridge- type 
IJJs, suggesting that ATh for sample S302 is ^0.23 K for 
Pdis =0.5 /iW in our measurements at the lowest temper- 
ature. These estimations also support that our results 
cannot be explained only by the self-heating effects. 

In order to check the influence of the self-heating in a 
more quantitative manner, we measured P(/) for the sec- 
ond switch by varying the values of dl/dt and the wait- 
ing times between the current-ramp pulses, as shown in 
Figs. [3a) and[S]Jb), respectively. Note that the changes 
in dl/dt correspond to those in the pulse widths since 
the amplitude of the ramped current pulse is fixed in the 
present measurement system. We measured the bridge- 
type IJJs (sample S302) for this purpose since the ef- 
fect of the self-heating process was expected to be more 
clearly observable. The experimental results indicated 
that 1/t(/) for the second switch was independent of the 
pulse widths and the waiting times. In terms of self- 
heating, this result can be understood if the experimen- 
tal time scale is longer than a time scale of the crossover 
from the short-time limit to the long-time limit since the 
temperature rise is already saturated at the current of 
the second switch.— However, even in such cases, Tosc 
should depend on the bias current, /, which determines 
the power density of the self-heating, giving rise to a 
convex curve in log-plots of t~^{I). In other words, Tesc 
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FIG. 7: (Color online). -P(/) for the second switch of sam- 
ple S302 under microwave irradiation at several frequencies, 
which was measured at 0.42 K. 

should be roughly proportional to or at least to a 
power of / larger than one. As shown in Fig. El^a), the 
measured data was compared with the calculated t^^(/), 
where Tosc was assumed to be constant (solid line), pro- 
portional to / (dotted hue), or proportional to P (dashed 
line), while was fixed to the value listed in Table I. 
We found that it was possible to successfully fit t~^(/) 
with the solid line rather than the dotted or the dashed 
lines over the whole range of the current distribution 
measurements. These results strongly suggest that the 
saturated behavior of Tosc below T^^""^ must be explained 
by considering the contribution of effects other than the 
self-heating effects in the resistive state after the first 
switch. 

Perhaps, this unusual behavior implies an increase of 
the electron temperature due to the injection of quasi- 
particles in the resistive state after the first switch, as 
previously discussed in terms of the modified heating the- 
ory of non-equilibrium superconductors/^ In this theory, 
the injected quasiparticles are assumed to interact with 
phonons with an energy greater than 2A, where A is the 
superconducting energy gap. However, this mechanism 
seems to be unable to explain the increase in the elec- 
tron temperature, which is as large as '^lO % of Tc, as 
discussed in the Appendix. Furthermore, in the mesa- 
type IJJs with a few junctions, the second switch should 
occur at a junction close to that at which the first switch 
occurs. On the other hand, in the case of the bridge-type 
IJJs, the junction at which the second switch occurs is 
not necessarily close to the junction at which the first 
switch occurs. Therefore, the mechanism of quasipar- 
ticle injection is expected to be different between both 
types of IJJs. Our results contradict this expectation. 

Another possible candidate is the MQT process for the 
second switch. In this scenario, the disagreement be- 
tween T^^nd g^j-^^ T'^f^ should be resolved self-consistently. 
Note that r^/"^ is estimated on the basis of the single 
junction model described in Sec. [Tll where the interlayer 



coupling effect is ignored. Thus, a more general treat- 
ment of Tcr in the multi-junction system is essential for 
resolving this issue. 

In order to examine the possibility for the occurrence of 
MQT for the second switch, wc also attempted to observe 
a resonant peak in P{1) for sample S302 by microwave 
irradiation. As shown in Fig. [7l wc did not observe any 
resonant peaks in P(/), although the distribution of P(/) 
was slightly shifted to the lower current region due to the 
strong microwave irradiation. However, these results do 
not rule out the possibility of MQT for the second switch. 
If the crossover temperature for the second switch is given 
by T^^"^, the energy- level spacing Eqx seems to exceed 
1 THz, suggesting that an n-photon process (n > 20) is 
required in order to observe ELQ. It is almost impossible 
to observe such higher-ordered multi-photon processes in 
the present experimental setup. Rather, the possibility of 
ELQ for the second switch should be tested by applying 
direct THz radiation, which will be performed in the near 
future. 



V. CONCLUSION 

We studied the MQT behavior of IJJs in BSCCO sin- 
gle crystals by comparing a small mesa, which is a few 
nanometers thick and consists of two or three IJJs, with 
narrow bridge structure consisting of a few hundreds of 
IJJs. We carefully distinguished the switched junction 
from the remaining junctions in the repeated measure- 
ments in order to obtain the switching current distribu- 
tion for both types of IJJs. For the first switching from 
the zero- voltage state, we observed the following three 
distinct regimes in the temperature dependence of P{I) 
for both types of IJJs: (i) the phase-retrapping regime 
at high temperatures, (ii) the TA regime at intermedi- 
ate temperatures, and (iii) the MQT regime at low tem- 
peratures. In the TA regime, the measured P{I) for the 
mesa-type IJJs were in excellent agreement with the con- 
ventional theory of a single Josephson junction. Together 
with the observation of a multiphoton absorption process 
between the quantized energy levels at the lowest tem- 
perature, we confirmed that the observed MQT behaviors 
were intrinsic to a single IJJ in BSCCO and independent 
of the device structure. 

For the second switch (from the first resistive state), 
we found that the temperature-independent behavior of 
Tesc, which is similar to the MQT behavior in the first 
switch, survived up to a much higher temperature than 
the crossover temperature estimated from the conven- 
tional single Josephson junction model. Based on our 
comparative study of two types of IJJs, we concluded 
that this unusual behavior for the second switch could 
not be explained in terms of the self-heating in the resis- 
tive state after the first switch. As possible candidates for 
the explanation of these observations, the MQT process 
for the second switch and the effective increase of elec- 
tronic temperature due to quasiparticle injection were 
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discussed. Unfortunately, the origin of this unusual be- 
havior remains unknown at present. Although the MQT 
behavior for the first switch of IJJs in the BSCCO sys- 
tem is well established in this work, the understanding of 
the second switching phenomenon remains incomplete. 
The following two approaches are required in order to 
obtain an important key to resolving this issue. One is 
the full treatment of the MQT rate considering the in- 
terlayer coupling effects in the stacked IJJs. The other is 
the direct observation of a resonant peak in P{I) under 
THz irradiation. The multi-dimensional aspects and the 
interlayer coupling effects are prominent features of IJJs 
systems in high- Tc cuprates. We believe that the phase 
dynamics of this system will attract considerable inter- 
est in the exploration of new physics in multi-junction 
systems and the application of phase qubits to future 
quantum information processing. 



as follows, 



' N 
N, 





0^ 











Nt ' 



2N^ 



-lo, 



(A.3) 
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where Nt{— ^/2NZt/Rt^) is the number of thermal 
equilibrium QPs, and Tfi{— 1/RNt) is the characteristic 
time at which QPs recombine into Cooper pairs. 

Eqs. ([Ol) and (AA\i suggest that NJN^t ^ {N/Nxf 
in the limit that r-y 3> t^, where the pair-breaking by 
phonons was more dominant than the other relaxation 
process of phonons. In the T*-model by Parker, the ef- 
fective temperature, T* , which characterizes the distribu- 
tion function of QPs through that of phonons, is related 
to the ratio {N/Nt), 
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APPENDIX 

In the modified heating theory proposed by Parker, 
it is assumed that the injected quasiparticles (QPs) in- 
teract with phonons with energy greater than 2 A, where 
A is the superconducting energy gap. The energy dis- 
tribution of non-equilibrium QPs is determined by that 
of phonons, which is described by a thermal distribution 
with an effective temperature T*. Thus, the properties 
of non-equilibrium QPs and phonons are described by a 
set of rate equations proposed by Rothwarf and Taylor, '^^ 



dN 
Itt 

dt 



2N^ 

TB 



2 TB ''7 



(A.l) 
(A.2) 



where N is the number of QPs, is the number of 
phonons with an energy greater than 2A, Iq is the rate 
of externally injected QPs, i? is a recombination coef- 
ficient, tb is the characteristic time at which phonons 
create QPs by a pair-breaking process, and is the char- 
acteristic time at which phonons disappear by processes 
other than the creation of QP. N^t is the number of ther- 
mal equilibrium phonons with an energy greater than 2A. 
The steady-state solutions of the rate equations are given 



\NtJ 
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(A.5) 

where X (or X*) is 2A{T)/kBT (or 2A{T*)/kBT*). 
Since the integral part included in Eq. (jA.sp is less 
than unity, {N/Nt}'^ gives the upper limit of the ratio 
{T*lTf. In the Umit of Tr/rB > 1, the following rela- 
tion is suggested. 



< 1 



laTR 



TB Nt 



= 1 



tbN^R 



(A.6) 



Here, we used the relation tr — 1/RNt- The inverse of 
the recombination coefficient, 1/i?, and the phonon pair- 
breaking time, Tb, for d-wave superconductors were mi- 
croscopically estimated by J. Unterhinninghofen et al^ 
The estimated values of 1/i? and tb for rf-wave supercon- 
ductors were ^10 ps and ^1 fs, respectively. A naive es- 
timation of using the thickness of the superconducting 
layer (~ 0.3 nm) and the velocity of sound (~ 10"^ m/s) — 
suggests that 1 x 10~^^ sec. We also estimated the 
number of thermal equilibrium QP's, Nt, as Nt ~ 3x 10'' 
by using the relation Nt ~ aim* /e^T in the microwave 
conductivity, CTi - 5 x 10*^ f^-^m'^ [lO], 1/t- lxlO"per 
seo^, and the junction size ^ 1 /imxl /^mxl.5 nm. In 
our experiments, if we take the values of /o(= I^^'^/e) ^ 
1 X 10'^'' per sec and T=0.4 K, the effective tempera- 
ture is required to be T* < 0.4 -I- 1 x 10~^ K which is 
much smaller than the T^^"''. Thus, it seems that our 
experimental result of high T^^'"^ cannot be explained by 
the increase of the electron temperature due to the injec- 
tion of quasiparticles in the resistive state after the first 
switch. 

So far, we have discussed the pair-breaking process by 
phonons. It is argued that the pairing mechanism of high- 
Tc cuprates might be derived from the antiferromagnetic 
spin fluctuation rather than the electron-phonon inter- 
action. Even in that case, essential part of the above 
discussion is valid, independent of the mechanisms of the 
pair-breaking. 
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